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Abstract 

We study waves in a chain of dispersively coupled phase oscillators. Two approaches - a quasi- 
continuous approximation and an iterative numerical solution of the lattice equation - allow us to 
characterize different types of traveling waves: compactons, kovatons, solitary waves with exponen- 
tial tails as well as a novel type of semi-compact waves that are compact from one side. Stability 
of these waves is studied using numerical simulations of the initial value problem. 
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I. INTRODUCTION 



The topic of this paper unifies two principal directions of nonlinear science: 
coupled self-sustained oscillators and soliton theory. Coupled autonomous self- 
sustained oscillators appear in different fields of science, they demonstrate a vari- 
ety of fascinating phenomena. In this study we demonstrate that particularities 
of the coupling for a rather simple setup - all oscillators are identical, periodic, 
and form a regular lattice with a nearest neighbor coupling - can lead to highly 
nontrivial wave structures. Remarkably, a dispersive coupling of the phases of 
dissipative oscillators results in a simple lattice equation which is equivalent to 
a Hamiltonian one. We show that different types of traveling waves exist in 
this lattice: compactons (ultra- localized waves with a compact support), usual 
solitary waves with exponential tails, as well as the corresponding kink-type 
solutions. These waves appear from rather general initial conditions and exist 
for long times. After very long transients they are destroyed due to inelastic 
collisions and evolve into phase chaos. 

Coupled autonomous oscillators are subject of high interest in nonlinear science [2]. 
When the coupling of the oscillators is weak, they can be described in the phase approx- 
imation [3j, where only a variation of the oscillator phases matters. The corresponding 
models are used for the description of lattices [U El H3 IE], globally coupled ensembles and 
networks 0, El EDI El 12] ■ In the absence of coupling, the phase equations have only zero 
Lyapunov exponents, therefore whether the phase dynamics is dissipative or conservative 
depends solely on the properties of the coupling. In studies which focus on synchronization 
properties, one typically assumes that the coupling is dissipative which thus tends to equalize 
the phases. However, certain types of coupling lead to a conservative dynamics - a prominent 
example is a splay state in a globally coupled ensemble of oscillators [131 EH E51 EH1 HTj . 

In the present work we consider the dynamics of a one-dimensional lattice of oscillators 
with a dispersive, conservative coupling. A realization of such a lattice may be a multi-core 
fiber laser [T5] , where individual self-oscillating lasers are arranged in a ring, or an array 
of Josephson junctions. Since both, the local phase dynamics and the coupling are non- 
dissipative, the dynamics is expected to be similar to that of the well-known Hamiltonian 
type. An example of a Hamiltonian lattice is the sine-Gordon lattice, for which the basic 
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building blocks are traveling solitary waves like pulses or kinks, that on integrable lattices 
collide elastically, and in the non-integrable cases evolve into chaos [19]. 

In recent years two new concepts appeared that have significantly extended our under- 
standing of nonlinear regimes in Hamiltonian lattices. One concept introduces localized 
periodic breathers [20]; the other introduces compact excitations in genuinely nonlinear lat- 
tices and wave equations. Unlike the usual solitons that have exponential, or algebraic, tails, 
the corresponding traveling waves have compact or almost compact support. These waves, 
named compactons, have been introduced in [2"T1 12"2l |2"3] . 

In the present paper we study general traveling waves in a chain of dispersively coupled 
nonlinear self-sustained oscillators, continuing our previous works (2U E5] focused on com- 
pacton solutions in phase oscillator chains. Contrary to [231 EH] , we consider here a generic 
phase model without additional symmetries. We show that dispersively coupled oscillators 
possess not only compactons, but also "classical" solitary waves with exponential tails. In 
dependence on parameters of the coupling and on the velocity these solution bifurcate. Re- 
markably, some of the compact solutions described in [2U ES] exist only in the symmetric 
situation, while other can be continued to the non-symmetric case. In the general, non- 
symmetric case, we describe a novel type of semi-compact waves which are compact from 
only one side and have an exponential tail on the other side. Another novel feature appears 
in the analysis of traveling waves on the lattice: here we describe a particular bifurcation 
from monotonic to oscillatory tails of a solitary wave, such a bifurcation does not exist in the 
quasicontinuous approximation where waves are described by means of a partial differential 
equation. We also report on a remarkable property of the solitary waves with exponential 
tails: they are extremely stable to collisions, so that a chaotization in a finite lattice occurs 
- contrary to the case of compactons - after exponentially long transients. 



II. BASIC MODEL 

We consider the phase dynamics of a chain of identical self-sustained oscillators with 
frequency uj 

^ = u + q((p n -i - V?n) + <?On+l - (p n )- (!) 

Here ip n is the phase of the n-th oscillator and q(ip) = q(<p + 2n) is the coupling function. 
Hereafter, we assume that this function is even, what corresponds to dispersive coupling, 
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see [25j for details. Introducing the phase differences v n 
rewritten as 

" V rf g(w) = q{v n+ i) - q(v n -i). 



dt 



ip n equation Q can be 



(2) 



The simplest possible even 27r-periodic coupling function is q(v) = cosv; this case was 
studied in [25], where it was shown, that compactons and kovatons (glued together compact 
kink-antikink pairs) are solutions arising out of the background v* = 0. Here, we extend 
the results of [25] and demonstrate that this lattice also bears periodic waves and solitons, 
as well as semi-compact waves. In [2S] it was also shown, that Eq. (JIJ is a Hamiltonian 
system. It possesses several integrals but is non-integrable. Additional symmetries in q(v) 
may result in additional symmetries in Eq. ([2]), we will encounter them below. 

Our main analytic tool for the study of the wave solutions of ([2]) is the quasi-continuous 
approximation (QCA). Associating n with a continuous variable x, one develops q{v n ±\) by 
the Taylor series expansion up to the third order 



Q{Vn±l / 



d 1 d 2 1 d 3 
dx 2 dx 2 6 dx 3 



q{v) 



(3) 



Inserting this into ([2]) results in the partial differential equation 

1 d 3 



dv 
di 



d_ 

dx 



6 dx 3 



q{v) 



(4) 



which is the QCA for the lattice. We stress here, that this approximation is not based on 
a small parameter: because spacing between the sites is one, the higher-order terms are in 
general of the same order as the lower-order ones. Thus the validity of this approximation 



may be supported only by a comparison of its predictions (Section III) with the numerical 
solutions of the full equations (J2l) (Section IV). 



III. TRAVELING WAVES IN THE QUASI-CONTINUUM 

In this section we analyze waves in the QCA Q. First, we note that any constant v = v * 
is a solution. We will look for waves on the base of such a flat profile. Inserting the traveling 
wave ansatz v(x,t) = v(x — \t) = v(s), where A is the wave velocity, into Eq. Q and 
integrating once yields 

X(v - v*) + 2(q(v) - q(v*)) + ~q(v) = 0. (5) 
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The integration starts at s , where it is assumed, that v(sq) = v * = const. This is not 
valid for periodic waves, where one has to introduce the curvature of q(v (sq)); however the 
resulting equations are equivalent to those derived below. We multiply Eq. ^ with dq(v)/ds 
and integrate again, to obtain 



A 



q(v) (v - v*) - Q(v, v*) + (q(v) - q(v*)) 2 + 
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dq 
ds~ 



1 2 



0. 



The function Q(v,v*) is defined as 



Q(v, u) = q(x)dx 



Equation (JoJ) can be rewritten as 



dq 
dv 



1 f dV Y TT ( S 
+U{V) 



with the potential 



U(v) 



X[q(v)(v - v*) - Q(v,v*)} + (q(v) - q(v*)f 



(6) 



(7) 



(8) 



(9) 



(q'(v)y 

Besides Eq. (|sj) one can also derive a system of first-order ODEs from Eq. ([5j, which reads 

dv du 3\(v-v*) + 6(q(v)-q(v*))+q"(v)u 2 , 

ds ds q'yv) 

In the following we also need the properties of the linear approximation q(v) = qiy*) + 



q'{v*){v — v*), then the equations <\W) simplify to 

du _ 3 A + 6q'(v*) 
ds" ~ ~ 



dv 



u 



(v-v*), 



(11) 



ds ds q'(v*) 

provided that q'(v*) 7^ 0. The stability of the fixed point at v — v* is determined by the 
eigenvalues of the Jacobian 



h,2 = ±v / -{3\ + 6q'{v*))/q'{v*). 



(12) 



A. Solitary waves 

Solitary waves on the base of the constant field v = v* are the homoclinic orbits of Eq. <\8 



and Eq. (10). They start at v*, grow to a peak at v m , and then go back to their origin v*. 
An equation for the wave velocity As can be obtained from the condition U(v m ) = which 
immediately yields 

A 5 = (g(^) -g(^*)) 2 / 13 x 
5 Q(v m ,v*) - q(v m )(v m -v*)' 
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Solitary waves with exponential tails. For the existence of a homoclinic trajectory 



one needs the fixed point of Eq. (10) to have one stable and one unstable direction, hence 
v* has to fulfill 

- 3A+ f^'-> > 0. (14) 
q'(v*) y ' 

This condition yields a critical velocity Ac = —2q'(v), which separates a saddle-type sta- 
tionary solution from a center. 

A solitary wave with exponential tails is shown in Fig. [4] below. The coupling function 
is q(v) = cos(v) and the background is v* = tt/4. The wave velocity is a free parameter, 
but bounded by the above condition. For the case v* = tc/A and v m > v* this results in 
A > 2sin7r/4 = \/2. In Fig. [Z]the wave velocity was chosen to A = tt/2. The tails of the 
solitary wave decay exponentially, corresponding to the eigenvalues of the stationary state 
v*. 

Solitary waves with compact tails - Compactons. Compactification may occur, 
if q'{v*) = 0. Linear waves do not exist around this point and the approximation in Eq. (11) 
does not hold. Therefore, one needs to approximate q{v ) to the second order in the vicinity 
of v*. When this approximation is inserted into Eq. ^ one obtains the K(2,2) equation, 
defined in [21]. For q(v) « q(v*) + av 2 a solution of (pj is given by 

{-?r cos 2 ( A /|s ) |x|<7rw| 
3a VV § J 1 1 V 3 (15) 

else. 

Usually, one cannot match two different solutions of one ODE, but here, the highest order 
operator degenerates at v — v* and the solutions uniqueness is lost. In the surrounding off* 
the solution will behave like the K(2, 2) and compact waves occur in the full phase equation 
Q. In Fig. [5] below we show a compacton for the coupling function q(v ) = cost> and v* = 0. 
The wave velocity was chosen to A = 2/n. 

B. Kinks 

The second class of traveling wave solutions are kinks, which correspond to heteroclinic 
orbits between two uniform states v* and v*. The height v* of the kink has to fulfill condition 



(13) and furthermore it has to be a fixed point of @, meaning that U'(v*) = 0. This gives 



the following condition for the speed of the kink (it can also be derived from ([5]) where one 




FIG. 1: The values of uniform background states that can be connected by kinks, in the QCA for 
the particular coupling q{v) = cosv + acos2v and various values of a. 

assumes a constant solution v m = const): 

x K = 2 m=m. (i6) 

v — V 



Combining the velocity relations (13) and (16) by As = Xk results in the final condition for 
the height of the kink 

Q(r,^) = wr)+ " ( ;' ))(r -" ,) . (it) 

In Fig.[T]the values of possible uniform states connected by a kink are shown for the particular 
coupling function q(v) = cosv + acos2t>. Note, that a bifurcation occurs at a = ±1/4 and 
two new branches of kink points emerge. For q(v) = cosv one finds v* = it — v* and the 
velocity belonging to this kink is X max = X K = Acosv* /(n — 2v*). 

Kinks with exponential tails. For kinks with exponential tails the same condition 



(|14j) as for the solitary waves has to be fulfilled. In Fig. |6j a kink with exponential tails is 
shown. The coupling function is q(v) = cos(w) and v* = ir/4, v* = 3tt/4. The velocity is 
A = V32/tt. 

Compact kinks - Kovatons. If v* and v* fulfill the compactification condition 
q'(v*) = q'(v*) = 0, both tails will become compact and may form a compact kink-antikink 
pair, named kovaton [21]. An example of this wave form is shown in Fig. [7j with q(v) = cosv 
and v* = 0, v* = tt and X K = 4/-7T. 

Exponential-compact kinks. In addition to kinks with exponential tails and kinks 
with compacton tails one can also observe semicompact kinks with one exponential and one 
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compact tail. Consider the coupling q{v) = cosv + acos2t> with a = 0.2. In this special setup 
v* — fulfills the compactification condition and v* = 2.39955 is the kink point satisfying 



(|17j) with velocity (|16j) \ K = 1.60011. This kink is shown in Fig. |8j It is compact at v = v 
and exponential at v = v*. 



C. Periodic waves 



Periodic waves around v* exist if the eigenvalues in (12) are purely imaginary. A 



periodic wave in the QCA is shown in Fig. [9j The velocity of a periodic wave must satisfy 

3A + 6q'(v* 



the condition resulting from (10): 



q'iy*) 



< 0. 



(18) 



For q(v) = cosf and < v* < n this condition yields A < 2sint>*. At A = Ac = —2q'(v*) 
some kind of bifurcation occurs. 



To quantify the dynamical behavior of the QCA near Ac we simplify (10) to 

dv 



ds 



u 



du 

d7 



X(V — V ) — V — V . 



(19) 



This ODE is not an approximation of (10) in a strict sense, since we have neglected the 



terms q"{v)u 2 and q'(v), but the qualitative behavior does not change. System (19) has two 
fixed points (vi,ui) = (v*,0) and (^2,^2) = (X — v*,0). The eigenvalues of the corresponding 
Jacobian are 



45 = ± V\ - 2v* if I = ±V2v* - A. 

If the velocity reaches the critical value Ac = 2v*, the two fixed points coincide and further- 
more both fixed points change their roles; if A < Ac the center is (vi,Ui) and (^2,^2) is a 
saddle-point; if A > Ac the roles have changed and (vi,ui) is the saddle and (t>2,«2) is the 
center. Now we return to Eq. (10). In Fig. [2] the phase space of (10) near the critical velocity 
is shown. If A < Ac the first fixed point will be a center, hence periodic waves exists in the 
QCA and for A > Ac the fixed point is a saddle point. Furthermore, the unstable and the 
stable manifold of this saddle are connected - a homoclinic orbit exists, resulting in solitonic 
solutions in the QCA. So, at the critical velocity we have a transition or a bifurcation from 
periodic waves to solitons. 
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FIG. 2: Phase space of ([10]) for v* = vr/4 and a) A = X c - 0.1 and b) A = X c + 0.1. The critical 
velocity is Ac = 2sin7r/4 = a/2. It is clearly visible that below the critical velocity the fixed point 
(t>i,ui) is a center, hence periodic waves exist in Q. Above Ac* one can observe solitonic solutions: 
the background v\ is a saddle point and a homoclinic orbit exists. 



IV. TRAVELING WAVES IN THE LATTICE 

Now we turn our attention to traveling waves in the full phase lattice model 1J2J. The 
wave ansatz for this model can be formulated as 



v n (t) = v(n — At) = v(s) 
with velocity A. Inserting (20) into ([2]) yields 

v = j(q{v(s-l)) -q{v(s + 1))). 
We integrate this equation from sq to s to obtain 



v (s) - v* = - J [q(v*) - q(v(r))]dn 



(20) 



(21) 



(22) 



where it is supposed, that v(s < sq) = v*. Again, as in the continuous version, the exact 
initial conditions v(so) are not relevant, they can be absorbed into the constant of integration. 
If one requires that v(s) = v* is a constant solution, hence a kink exist, the corresponding 
velocity X K has to satisfy 



A 



K 



.q(v*)-q(v*) 



v — V 



(23) 



This condition is exactly analogous to condition (16) in the QCA. 
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A. Fixed point analysis 



Similar to the fixed point analysis in the QCA, one can analyze the behavior of traveling 



waves close to the background v*. To this end we linearize q(v) « q(v*) + a(v — v*) in (21) 
and apply the exponential ansatz v(t) = Aexplt. This yields the characteristic equation 

l = \{e- l -e l ). (24) 

Note again, that a = q'(v*) ^ 0, meaning that this approximation does not hold for the 
compacton backgrounds. 

We split I into its real and imaginary part I = p + iq to obtain 

p = — 2 — cosgsinhp and q = 2— sing cosh p. (25) 
A A 

For a purely imaginary eigenvalues p = we obtain 

a sing 

g = — 2— sing or A = — 2a . (26) 

A g 

This function is plotted in Fig.|3](a). In this plot, the dots mark possible points for transitions 
to eigenvalues with real parts. In Fig. [3](b) all eigenvalues I = p + iq are shown. Purely real 
eigenvalues are 

a . „ sinhp . . 

p=— 2— smhp or A = — 2a . (27) 

A p 

So, when A crosses —2a (point 1 in Fig. |3](a)), a bifurcation from two purely imaginary 
eigenvalues to two purely real eigenvalues occurs. This scenario corresponds to the transition 
from periodic to solitary waves and the critical velocity is Ac = —2a. The situation is 
analogous to the bifurcation in the QCA and the critical velocity is the in both situations. 

The next bifurcation occurs, when A crosses point 2, see Fig. |3j Then, a bifurcation from 
a center to a stable and an unstable spiral point occurs. This refers to the transition from 
periodic waves to solitary waves with oscillating tails and exponentially decaying amplitude. 
Since the bifurcation occurs on the imaginary axis, one can calculate the critical velocity from 



(25) by setting A'(g) = and for the special case q(v) = cosv one obtains A m —2a ■ 0.217. 
Note, that there is no counterpart in the QCA for this bifurcation. 
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(a) 



(b) 



FIG. 3: (a): Wave speed A in dependence of the imaginary eigenvalue q, see (26). The dots mark 
possible bifurcation points to real eigenvalues, (b): Wave speed A vs. the real and imaginary part 
of the eigenvalues, with a = sin(0.2). Only the positive parts of the real and imaginary axes are 
shown. 

B. Numerical determination of traveling waves 



From (22) it is possible to construct a numerical scheme to find solitary wave solutions 



of the lattice [251 1211 [27] . One initially guesses a wave profile v (t) and than iterates 

1 r s+1 f \\v 1 1 \ 3 ^ 2 

v(s) = v * + xJ x ^ ~ <?W r ))) dr = (w ) ^ (28) 

|| ■ || denotes the Li-norm. The integral is calculated by a high order Lagrangian integration 



rule [28j. To construct kink solutions, one has to omit the normalization in (28) by setting 



Vk+i = v. Periodic waves can be obtained by a slight modification of (28). Here, the wave 



length w is introduced and periodic boundary conditions v(0) = v(w) are assumed in (28). 



We want to point out two issues one has to keep in mind when using this algorithm. 



First, in (28) the normalization exponent 3/2 is used. In a few cases this exponent is to 
large and has to be set to smaller values, otherwise the algorithm will diverge. Secondly, for 
backgrounds different from one has to shift the coordinates v > v* + v. 

C. Solitary waves 

Solitary waves arise out of a background with one stable and one unstable direction. So, 
they have to fulfill (A + 2q'(v*))/q'(v*) > in order to obtain two real eigenvalues of the 
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FIG. 4: Top panel: Soliton solution for the coupling q{y) = cos v, which arise out of the background 



v* = 7r/4. The form of the soliton was calculated numerically with the help of (28 ) and the velocity 
is A = 7r/2. Here and in the following figures bold dots show the soliton on the lattice and the 
dashed line the solution of QCA ^ which has an additional offset for better visibility. Bottom 
panel: the soliton in logarithmic scale. 

fixed point. 

Solitary waves with exponential tails. In Fig. [4] we show a solitary wave with 
exponential tails. The coupling is q(v) = cosv and the background is v* = tt/4. The 
velocity of was chosen to A = 7r/2, fulfilling the condition (A + 2q'(v*))/q'(v*) > 0. The 



solitary wave was computed with the scheme (28) 



Compact solitary waves. Again, as in the QCA, the condition q'(v*) = allows the 
compactification of the tails. In Fig. [5] we show a compacton arising out of the background 
v* = for the coupling q{y) = cosv. The wave velocity was set here to A = 2/ir. The 
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FIG. 5: Top panel: the shape of the compacton for the coupling q{v) = cost;. The background is 
v* = and the velocity was set to A = 2/tt. Markers show the compacton on the discrete lattice 
and the dashed line is the corresponding solution in the quasi- continuum approximation. Bottom 
panel: the same plot in logarithmic scale. The super exponentially decaying tails are clearly visible. 

compacton is not purely compact, but has super-exponentially decaying tails [25j. Thus, 
although there is a qualitative difference between the lattice and the QCA, quantitatively 
these solutions are very close to each other. 



D. Kinks 



To derive an analogon to (17), we rewrite Eq. (23) as Ak(v* — v*) = 2(q(v*) — q(v*)). 
multiply with q'(v*), and integrate over v* from v* to v* to obtain 



Xk I q'(v ) (v — v*)dv = 2 / q'(v) yq{v*) — q(v) jdv 



(29) 
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and finally 



X K = (g(^)-g^*)) 2 (30 ) 
Q(v*,v*) — q(v*)(v* — f*) 



Combining Eq. (23) and (30) yields 



Q(r,^ )= fa(r)+ ^; ))(r -"' ) , pi) 

which matches exactly the kink condition for the QCA. 

Kinks with exponential tails. In Fig. [6] we show the shape of a kink. The coupling 



is q(v) = cosv and the background is v* = it/ 4. The velocity of the kink is given by (23) 
X K = V32/tt. 

Kinks with compact tails. One can observe compact kinks, if q'(v*) = and q'{v*) = 
0. For the coupling q(v) = cosv such a case exists with v* = and v* = it. The shape of 
this compact kink is shown in Fig. [7j Here, the velocity is A = A/tt. 

Semi-compact kinks. It also possible to observe kinks with one exponential decaying 
tail and one compact tail. This is the case for q{v) = cosv + acos2t> with a = 0.2. For 
A = 1.60011 and v* = 2.39955 the kink condition ( |3l| is satisfied and such a kink is found 
by the numerical method described above, see Fig. [8j 

E. Periodic waves 



Periodic waves can be calculated with (28) and periodic boundary conditions. An example 
is shown in Fig. [9j Here, the offset is v* = ir/4, the wave length is w = 5tt and the velocity 
is A = 7r/2. 

F. Solitary waves with periodically decaying tails 



From the fixed point analysis of the advanced-delayed equation (21) a bifurcation occurs 
at point 2 in Fig. [3|a). So, if the velocity A reaches the critical point X C} the fixed point 
changes its type from a center to a stable and an unstable focus. This corresponds to a 



solitary wave with oscillatory decaying tails. In Fig. [T0]we show an example of such a wave. 
The offset is v* = —0.2 and the wave velocity is A = 1.0. This behavior does not occur in 
the quasi continuum. 
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FIG. 6: Top panel: The shape of the kink for v* = 7r/4, the velocity is A = y/32/ir. Center panel: 
The kink in logarithmic scale. Bottom panel: The kink shown from its top in logarithmic scale. 
Markers show the kink on the lattice and the dashed line is the QCA with an additional offset for 
better visibility. 
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FIG. 7: Top panel: The wave form of the compact kink solution between v* = and v* = n for 
the coupling q(v) = cos v. Markers show the wave form on the discrete lattice and the dashed line 
represents the quasi- continuum approximation of the kink. Bottom panel: The sine of the kink in 
logarithmic scale. 

V. NUMERICAL SIMULATIONS OF THE INITIAL VALUE PROBLEM 

In this section we demonstrate, how the traveling waves described in the previous sections 
appear in the course of the evolution of the lattice. We will restrict ourselves to the simplest 
coupling term q(v) = cosv, while the background state will be general v* ^ 0. Furthermore, 
we will study the stability properties of the colliding waves. 
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FIG. 8: Top panel: The shape of a kink with one exponential and one compact tail. The coupling 
for this specific wave is q[v) = cos v + a cos 2v with a = 0.2. The position of the kink point and the 



wave velocity can be obtained from (31 1 and (23), their values are A = 1.60011 and v* = 2.39955. 
Middle panel: The kink in logarithmic scale. The compact tail of the kink is hear clearly visible. 
Bottom panel: The kink shown from its top in logarithmic scale, the exponential tail becomes 
visible. 
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FIG. 9: The shape of a periodic wave on the lattice. The offset is v* = 7r/4, the wave length is 
w = 5w and the velocity A = n/2. 



A. Evolution of an initial pulse 

First, we consider the evolution of an initial cos-pulse with the coupling q(v) = cos v. 
The initial condition is 

+ + cos (*=>)] \n-n \<w 
v n {0) = < (32) 

I v* else, 

where A is the amplitude, n is the center and w is the half width of the pulse. 

In Figs. [TTp2" we compare the evolution for different values of v*. In Fig. 11 we set v * — 



(q'(v*) = 0) and one can observe compactons and kovatons arising from the initial pulse. 
In Fig. |ITj(a) a wave train of compactons emerges out of the initial pulse. The speed of the 
compactons increases with increasing amplitude. In Fig. [TT](b) we have increased the width 



and the amplitude of the pulse and one kovaton is observed. In the bottom plot of Fig. [TT 
a narrow initial pulse creates a wave source, emitting periodic waves. 



In Fig. 12 we show results for the background v* = tt/4. Here, q'(v*) ^ and the solitary 
waves arising from the initial pulse possess exponential tails. Fig. [l2|a),(b) are similar to 
Fig |TT|(a),(b), where the initial pulse decomposes into a train of solitons and kink. Note, 
that the number of emitted solitary waves is smaller than for the case v* = 0. Furthermore, 
periodic waves around v* can emerge, see Fig. [l2|(c) and (d). The plot in (d) is somehow 



similar to Fig. jTlTc) , with the difference, that appearing periodic waves are around v*. 
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FIG. 10: Solitary wave with periodic and exponentially decaying tails, v* 



-0.2 and A = 1.0. The 



plot was generated with the help of the traveling wave scheme (28) for the lattice equation; it has 
no QCA counterpart. 



Fig. 12 e) shows the evolution of a narrow initial pulse with a relative large amplitude. It 



results in a kink with periodic waves around the top of the kink. A detailed analysis of all 
possible waveforms goes beyond the scope of this paper and will be reported elsewhere. 

B. Transition to chaos in a finite lattice 



Wave trains shown in Figs. Ilp2 are obtain for an effectively infinite lattice (during the 
calculation times the boundaries are not reached). In a finite lattice, collisions between 
waves occur. We have used periodic boundary conditions, and observed that at large times 



eventually a chaotic regime appears. In Fig. [13] we show the evolution of an initial cos pulse 
with v* = 0.1. The upper plot shows the initial decomposition of this pulse into one kink and 
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FIG. 11: Evolution of different initial pulses for the coupling q(v) = cosv and v* = 0. The initial 
conditions were set according to ( |32| ). (a) w = 10 and A = 2, (b) w = 15 and A = 3 and (c) ui = 1 
and A = 1, The dashed line shows the initial condition u n (i = 0) and the solid line the lattice at 
the time t = 100. 

several solitary waves. These structures appear to survive collisions quite unaffected. The 
lower plot shows, that after some transient time chaos emerges. The chaotic state begins 
to develop, after a collision of two solitons produces a large-amplitude soliton-antisoliton 
pair. Then an avalanche of soliton-antisoliton collisions is triggered on, resulting in a fast 
chaotization. 



In Fig. [14] we show a remarkable dependence on the average transient time, after which 
chaos establishes, on the parameter v*. For larger values of v* the transient time is exponen- 
tially large, what means extreme stability of the solitary waves. Qualitatively, this stability 
can be attributed to a smallness of effects of discreteness of the lattice for large v*. Here, the 
waves are relatively wide, thus they are well approximated in the QCA, which is close to the 
integrable Korteveg-de Vries equation. For small v* the waves are close to compactons that 
are short and for them the discreteness that causes non-elasticity of collisions is essential. 
Furthermore, the number of emitted waves decreases with increasing v* and the velocity of 
the waves is bounded from below. These two effects reduce the possibility that two waves 
meet each other, resulting in an increased transient time. 
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FIG. 12: Evolution of different initial pulses for the coupling q(v) = cosv and the background 
v* = 7r/4, The initial conditions were set according to ( |32| ). (a) w = 10 and A = 1, (b) w = 15 and 
A = 1.5, (c) u; = 10 and A = 2, (d) w = 1 and A = 1 and (e) w = 1 and A = 2. The dashed line 
shows the initial condition u n (t = 0) and the solid line the lattice at the time t = 100. Furthermore 
the position of the kink point at 3/47T is shown. 




FIG. 13: Transition to chaos. The background is v* = 0.1 and the lattice contains of N = 100 sites 
with periodic boundary conditions. The field is shown in a gray scale versus time (horizontal axis) 
and space (vertical axis). Upper plot: A kink and several solitons emerge out of an cos pulse. The 
time interval is < t < 400. Lower plot: Emergence of chaos after a collision of two solitons which 
creates an soliton-antisoliton pair. The time interval is 2600 < t < 3000. 
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FIG. 14: Transient times to chaos for different backgrounds v*. The length of the lattice is N = 32 
and the initial cos-pulse is v n = 7r/4(cos(7r(n — N/2)/w) + l) + v* for |n — N/2\ < w and v = v* 
else. In order to obtain an average of the transient times we also varied the width of the initial pulse 
from 5 to 15 and calculated the transient time as the average over the transient times for different 
initial pulses. The line is an exponential fit and the transient times scales with T ~ exp 16. 7v* . 

VI. CONCLUSIONS 

We have demonstrated a variety of nontrivial wave structures in dispersively coupled oscil- 
lator lattices. Remarkably, they appear in a very simple lattice described by Eq. ([2]). In this 
study we have focused, contrary to previous works [2H[25], on the features that appear for 
a general, non-symmetric coupling function. While some nontrivial solutions (compactons) 
survive in a general case, other (kovatons) exist only in the symmetric situation. Instead, 
for a general case we have reported a novel type of semi-compact waves. In our study of the 
waves on a lattice we have described a novel transition from monotonic to oscillatory tails of 
solitary waves that does not exist in the quasicontinuous approximation. Our comparison of 
general typical solutions of the lattice model with special ones studied in [2U [25] has shown 
that the waves with exponential tails are much more "resistant" to chaotization compared 
to the compactons. 

Here, we would outline several possible extensions of the analysis. In general, coupling 
between oscillators can possess both dispersive and dissipative parts. The waves described 
in this work will still be observed if the dissipation is sufficiently small, this is confirmed by 
the perturbation analysis in [25]. Another feature disturbing the waves is a non-homogeneity 
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of the lattice, e.g. due to non-uniformity of coupling. We expect that the waves will scatter 
on such inhomogeneities, but this issue has not been studied yet. Finally, it is intriguing, 
what kinds of waves can be observed in two- and three-dimensional lattices; results in this 
direction will be published elsewhere. 

We thank P. Rosenau for helpful discussions and DFG for support. 
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